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I.  INTRODUCTION 


For  many  electromagnetic  problems  that  require  a  rigorous  solution  for 
the  electric  field,  a  closed-form  analysis  is  impractical  because  of  the  com¬ 
plexity  of  the  integrals  involved.  On  the  other  hand,  the  computer  approxi¬ 
mation  of  these  problems  is  generally  easy  and  accurate. 

An  approximation  method  that  is  extensively  used  in  the  community  is  the 
method  of  moments.^  The  method  of  moments  reduces  a  nonlinear  equation  to  a 
set  of  linear  equations  that  can  be  solved  by  matrix  methods  such  as  Gaussian 
elimination  and  conjugate  gradients. ^ in  this  report,  the  method  of  moments 
is  applied  to  Pocklington* s  integral  equation  to  approximate  the  electric 
field  radiated  by  a  wire  dipole  antenna  two  wavelengths  long. 


II.  PROBLEM  FORMULATION 


In  order  to  compute  the  far-field  pattern  from  the  dipole,  the  procedure 
is  first  to  compute  the  current  on  the  wire  by  the  method  of  moments  and  then 
to  integrate  the  current  using  the  Green's  function.  The  method  of  moments 
can  be  used  to  approximate  the  current  on  the  wire  as  pulses  at  discrete 
positions  on  the  wire.  The  distances  between  these  pulses  are  arbitrary  but, 
for  programming  ease,  they  should  be  uniform.  The  model  that  is  developed  for 
the  dipole  antenna  is  a  wire  divided  into  an  Integral  number  of  segments  with 
half-segments  at  the  ends  where  the  current  is  zero.  The  pulse  approximation 
dictates  that  the  values  for  the  current  and  the  voltage  on  each  segment  be 
constant.  The  current  and  voltage  can  be  represented  as  twin  arrays,  with 
each  element  of  an  array  corresponding  to  the  value  either  of  the  current  or 
the  voltage  on  a  segment.  For  the  center-fed  dipole  antenna,  as  shown  in 
Figure  1 ,  the  values  for  the  elements  of  the  voltage  array  are  zero,  except 
that  the  element  corresponding  to  the  center  segment  has  the  value  of  the  vol¬ 
tage  applied,  VA.  The  values  of  the  elements  of  the  current  array  are  the 
unknowns  in  the  dipole  antenna  problem. 

•  v 

The  number  of  segments  is  critical  to  the  resolution  of  the  current.  A 
rule  of  thumb  is  to  select  at  least  10  segments  per  wavelength.  For  example, 
if  the  antenna  is  O.SX  long,  then  at  least  five  segments  are  required.  The 
overall  size  of  the  dipole  antenna  to  be  analyzed  is  limited  by  the  number  of 
segments  that  can  be  treated  practically.  For  antenna  structures  that  are 
more  general  than  a  simple  linear  dipole,  the  two-dimensional  current 
distribution,  rather  than  axial  dipole  current,  must  be  obtained.  Thus,  the 
overall  size  of  the  antenna  structure  is  limited  to  designs  that  are  small  in 
terms  of  wavelengths. 


Figure  1.  Center-Fed  Dipole  Antenna  with  N  Segments.  Array  elements 
represent  the  current  and  voltage  on  the  segments. 


10 


III.  METHOD  OF  MOMENTS 


The  pulse  model  of  the  dipole  antenna  is  Implemented  In  the  method  of 
moments  by  using  the  pulse  function 

11 

0 

as  a  basis  function.  Because  zn  denotes  the  center  of  the  nth  segment  and  A1 
is  the  length  of  the  segment,  the  pulse  function  is  zero  except  at  the  nth 
segment.  Other  basis  functions  that  are  appropriate  for  modeling  the  dipole 
antenna  are  the  sinusoidal  and  triangular  functions.  The  equation  to  be 
linearized  in  the  wire  antenna  problem  is  Pocklington’ s  Integral: 

2  L 

(it2  +  J  c(z')G<z,z")dz'  -  “:Ji<*weoE^an(z) 

Sz  o 

where  z' ,  the  integration  point,  is  the  position  of  the  current  on  the  wire; 
z,  the  observation  point,  is  the  position  of  the  impressed  electric  field  on 
the  wire;  and  G(z,z")  is  the  Green's  function  for  a  thin  wire  antenna  of 
radius  a: 


elsewhere 


G(z,z ")  -  exp{-jk_[a2  +  (z-z_-)2]1/2  } 

/a2  +  (z  -  z')2 

The  current  on  the  wire,  C(z"),  is  the  unknown  for  which  the  moment  method 
will  approximate  a  solution.  The  observation  points  are  selected  to  be  on  the 
surface  of  the  wire,  where  the  expression  for  the  Impressed  electric  field  is 
valid,  and  the  integration  points  are  on  the  axis  of  the  wire.  The  first  step 
of  the  moment  method  is  to  write  C(z')  in  terms  of  the  basis  function: 


II 


C(z') 


I  Cn  P<z-,zn) 
n-1 


If  we  bring  the  second  derivative  inside  the  integral  in  accordance  with 
Pearson,^  the  Pocklington* s  integral  becomes 


N 


I 

n-1 


C 


n 


z  +  A1 
n 


(k2G  +  ^-|)dz' 


-jo>4ire  E. 

o  tan 


This  simplified  equation  can  be  linearized  by  using  a  Dirac  delta  function  as 
the  testing  function.  This  changes  Pocklington* s  equation  to 


N 


I  c 


z  +  A1  , 
n  r,  2 


n-1 


2  G<VZ'>  +  [— 2  G<z*z')]2.2  \dz' 

z  -  Ll  52 

n  r 


-j  <d4 ire  E,.  „(z  ) 
o  tan  m 


The  second  derivative  can  be  evaluated  two  ways  -  analytically  and  by  finite 
differences.  The  difference  between  the  far-zone  electric  field  of  the 
analytical  and  finite  differences  is  minimal,  as  is  shown  in  the  numerical 
results  section  of  this  report.  However,  the  current  at  the  center  of  the 
wire  and,  therefore,  the  input  admittance,  do  not  agree. 

The  resultant  matrix  equation,  [Z][C]  -  [Al»E],  can  be  solved  for  the 
unknown  current  by  standard  matrix  inversion  methods,  such  as  Gaussian 
elimination  and  L-U  factorization  or  by  the  conjugate  gradient  method.  The 
conjugate  gradient  method  is  discussed  subsequently  because  it  is  not  as 
universally  used  as  the  other  methods  and  because  it  requires  less  memory. 
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IV.  CONJUGATE  GRADIENT  METHOD 


The  conjugate  gradient  method  is  an  iterative  method  that  assumes  a 
solution  and,  by  computing  the  residual,  adjusts  the  initial  solution  until 
the  relative  error  is  within  the  specified  tolerance.  The  method  is  best  used 
to  minimize  or  maximize  an  optimization  function  but  can  also  be  used,  as  in 
this  case,  to  solve  a  linear  equation.  The  conjugate  gradient  algorithm  is  as 
follows : Z 


zc0  -  V 


-b_,z*r 


Ck+1  *  °k+  tkpk 


rk+l  *  rk  +  ckZPk 


Pk+1  *  pk  "  \Z  rk+l 

tk«— 4 

,zpk« 


*  2 
,Z  rk+l' 


where 


V  -  [A1*E  ] 

Z*  ■  the  complex  conjugate  of  the  impedance  matrix 


The  algorithm  computes  the  residual  array,  rk,  at  each  iteration  and  makes 
adjustments  to  the  current  array,  C^,  based  on  the  residual  array  to  minimize 
the  relative  error.  The  expression  for  the  relative  error  is  lZCk  -  Vl/|v|. 
If  the  current  is  initially  assumed  to  be  zero  at  every  segment,  the  error 
equals  l  at  the  first  iteration  and  decreases  monotonically,  until  the  error 
is  less  than  the  prescribed  tolerance.  For  the  problem  of  the  wire  antenna 


two  wavelengths  long,  a  tolerance  of  10~°  was  used,  and  the  algorithm 
converged  within  a  reasonable  number  of  iterations.  A  chart  of  the  conver¬ 
gence  rate  is  shown  in  the  numerical  results  section  of  this  report. 

The  memory  requirement  for  the  conjugate  gradient  method  is  7N,  as 
opposed  to  2N^  +  2  for  the  IMSL  library  routine  that  was  used  to  verify  the 
results  presented  in  this  report. 


V.  NUMERICAL  RESULTS 


The  numerical  results  reported  in  the  tables  and  figures  of  this  section 
are  illustrations  of  the  comparison  of  four  combinations  of  second  derivative 
methods  with  linear  equation  methods:  (1)  analytical  derivative  with  con¬ 
jugate  gradients,  (2)  analytical  derivative  with  L-U  decomposition,  (3)  finite 
differences  with  conjugate  gradients,  and  (4)  finite  differences  with  L-U 
decomposition.  (The  listings  for  the  computer  programs  written  in-house  are 
given  in  the  Appendix.)  The  combinations  were  applied  to  a  theoretical  wire 
antenna  with  a  length  of  2X  and  a  radius  of  0.0I3X,  with  1  V  applied  to  the 
center.  Because  an  error  tolerance  of  10  was  set  for  the  conjugate  gradient 
method  and  10“^  for  the  L-U  decomposition,  the  current,  far-zone  electric 
field,  and  input  admittance  calculated  by  each  of  these  matrix  methods  for  the 
same  second  derivative  method  are  equal  to  at  least  eight  significant  fig¬ 
ures.  Because  the  numerical  values  calculated  by  the  conjugate  gradient 
method  and  L-U  decomposition  are  equivalent,  no  distinction  is  made  between 
them  in  Figures  2,  3,  and  4  or  in  Tables  1  and  2. 

Figures  2,  3,  and  4  illustrate  the  trends  of  the  current  approximated  by 
the  analytical  derivative  and  by  finite  differences  as  the  number  of  segments 
is  increased  from  10  to  20  and  then  to  60.  Assuming  that  the  current  approxi¬ 
mated  by  60  segments  is  close  to  an  actual  experimental  sampling  at  those 
segments,  we  can  deduce,  from  the  plots  of  the  approximations  with  10  and  20 
segments,  that  10  segments  is  not  sufficient  and  that  20  segments  is  suffi¬ 
cient.  In  fact,  if  we  were  able  to  superpose  the  plot  with  20  segments  over 
the  one  with  60  segments,  we  would  see  that  they  were  very  close  in  numerical 
value  but  that  the  one  with  20  segments  had  less  resolution. 

The  electric  field,  at  a  distance  of  50  wavelengths,  calculated  from  the 
approximations  of  the  currents  on  20  segments,  is  shown  in  Table  1.  The  dif¬ 
ferences  in  the  electric  fields  are  within  0.2  dB  for  the  main  lobes  and  2  dB 
for  the  low-power  angles. 
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Ten  Segment  Approximation.  Amplitude  and  phase  of  the  current  on 
a  center-fed  antenna  two  wavelengths  long,  0.U13  wavelengths  in 
radius  with  1  V  applied. 
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Twenty  Segment  Approximation.  Amplitude  and  phase  of  the  current 
on  a  center-fed  antenna  two  wavelengths  long,  0.013  wavelengths  in 
radius  with  1  V  applied. 
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Figure  4. 
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Sixty  Segment  Approximation.  Amplitude  and  phase  of  the  current 
on  a  center-fed  antenna  two  wavelengths  long,  0.013  wavelengths  in 
radius  with  1  V  applied. 


Table  1 .  Far-Zone  Electric  Field 


Angle 

(deg) 

Finite  Differences 
(dB) 

Analytical 

(dB) 

10 

-3.4 

-5.3 

20 

1.0 

+0.3 

30 

5.5 

6.3 

40 

11.7 

12.2 

50 

15.8 

15.9 

60 

17.0 

16.8 

70 

14.4 

14.2 

80 

5.7 

5.8 

90 

-2.3 

-3.8 

Although  the  agreement  of  the  far-zone  electric  field  is  adequate,  the 
admittance  calculated  with  analytical  derivatives  and  finite  differences 
yields  approximately  23  percent  difference  in  magnitude.  This  magnitude 
difference  for  the  2A  antenna  is  in  the  same  range  of  differences  reported  by 
Butler  and  Wilton^  for  a  1.5A  antenna.  The  value  of  the  input  admittance  is 
important  for  matching.  The  result  for  20  segments  are  given  in  Table  2. 


Table  2 .  Input  Admittance 


REAL 

Finite  differences 
Analytical 

IMAGINARY 

Finite  differences 
Analytical 


1.687  nmhos 
1.668  mmhos 


2  .466  inmhos 
1.609  mmhos 


The  calculations  were  executed  with  ANSI  Standard  Fortran  77  on  the  CDC 
Cyber  176.  Two  matrix  routines  were  used  to  find  the  solution  of  the  matrix 
equation:  (1)  the  conjugate  gradient  algorithm  described  in  this  report  and 
(2)  an  IMSL  library  routine  that  uses  L-U  decomposition  and  iterative 


improvement.  The  execution  time  varied  with  the  method  used  for  the  linear 
equation  solution  and  the  method  used  for  evaluating  the  second  derivative  of 
the  Green's  function.  For  the  second  derivative  evaluated  by  finite  dif¬ 
ferences  on  a  wire  approximated  by  20  segments,  the  times  for  different  linear 
equation  methods  are  given  in  Table  3 . 


Table  3.  CPU  Times  for  Finite  Differences 


Conjugate  gradient 
L-U  decomposition 


2.3  cpu  sec 
1 .9  cpu  sec 


The  times  for  different  second  derivative  methods  using  the  conjugate  gradient 
method  are  given  in  Table  4. 

Table  4.  CPU  Times  for  the  Conjugate  Gradient  Method 


Finite  difference 
Analytical  derivative 


2.3  cpu  sec 

5.3  cpu  sec 


The  times  given  in  Tables  3  and  4  are  the  total  cpu  times  for  calculating  the 
current,  input  admittance,  and  far-zone  electric  field.  Execution  times  for 
the  conjugate  gradient  method  were  made  with  an  initial  current  of  zero  at 
every  segment.  Because  the  conjugate  gradient  method  is  iterative,  the  exe¬ 
cution  time  depends  upon  the  residual  computed  from  the  initial  current  array. 

For  this  antenna,  for  which  the  current  is  symmetric  about  the  center 
segment,  the  conjugate  gradient  method  arrived  at  a  solution  in  N/2  itera¬ 
tions.  In  Figure  5,  the  logarithm  of  the  relative  error  is  shown  versus  the 
iteration  number  for  10,  20,  and  60  segments.  The  relative  error  approaches 
the  tolerance  of  10“®  monotonically  in  each  case. 
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VI.  CONCLUSIONS 


In  this  report,  the  numerical  methods  are  discussed  for  analyzing  a 
linear  antenna  using  the  method  of  moments.  Within  the  method  of  moments 
there  are  several  options  to  consider:  which  basis  and  testing  functions  to 
use,  how  to  evaluate  the  second  derivative,  and  which  linear  equation  method 
to  use.  The  criteria  used  in  selecting  the  options  may  be  as  follows: 

1.  Agreement  of  numerical  results  with  previously  accepted  results. 

2.  Ease  of  implementation. 

3.  Computer  execution  time. 

4.  Computer  memory  requirements. 

A  summary  of  the  options  discussed  in  this  paper  and  how  they  meet  the 
criteria  is  given  subsequently. 


A.  ANALYTICAL  DERIVATIVE  VERSUS  FINITE  DIFFERENCES 


The  far-zone  electric  fields  predicted  by  each  of  these  methods  for 
computing  the  second  derivative  are  in  relative  agreement,  but  the  input 
admittances  and  execution  times  differ.  The  finite  difference  calculation 
takes  less  cpu  time,  but  the  resultant  input  admittance  differs  by  23  percent 
from  the  admittance  calculated  with  the  analytical  second  derivative. 

B.  CONJUGATE  GRADIENT  VERSUS  IMSL  LIBRARY  ROUTINE 


The  results  of  the  conjugate  gradient  and  IMSL  routines  for  the  current, 
admittance,  and  far  field  are  equal,  but  the  routines  differ  in  memory 
requirements  and  execution  times.  The  IMSL  routine  takes  less  execution  time 
but  requires  more  memory  than  the  conjugate  gradient  routine. 


Numerical  results  in  this  report  should  be  valuable  in  matching  the 
options  to  specific  applications. 
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APPENDIX 


The  three  computer  programs  that  are  listed  compute  the  current 
distribution  and  far-zone  electric  field  for  a  center-fed  antenna  using 
Pocklington 's  equation  and  the  method  of  moments.  The  programs  run  inter¬ 
actively  with  voltage  source  amplitude,  phase  and  frequency,  antenna  length 
and  radius,  and  far-zone  distance  as  the  input  parameters. 

Listing  1  evaluates  the  second  derivative  analytically  and  solves  the 
matrix  equation  using  the  conjugate  gradient  algorithm. 

Listing  2  evaluates  the  second  derivative  by  finite  differences  and 
solves  the  matrix  equation  using  the  conjugate  gradient  algorithm. 

Listing  3  evaluates  the  second  derivative  by  finite  differences  and 
solves  the  matrix  equation  using  an  IMSL  routine. 


LISTING  NO.  1 


t  <;>u  c 

110  C  PROGRAM  WANTCG  CALCULATES  THE  CURRENTS  ON  A  WIRE  ANTENNA 

120  C 

130  C  WANTCG  CALCULATES  THE  CURRENTS  ON  A  WIRE  ANTENNA  WITH  THE 

140  C  SOURCE  AT  THE  CENTER.  (C)  IS  THE  CURRENT  VECTOR,  (V)  IS  THE 

150  C  APPLIED  VOLTAGE  VECTOR  AND  OPERA T  IS  THE  INTEGRAL 

160  C  OPERATOR  SUBROUTINE.  WANTCG  SOLVES  FOR  THE 

170  C  CURRENTS  USING  THE  CONJUGATE  GRADIENT  METHOD. 

180  PROGRAM  WANT CG < INPU T , TAPE5= I NPUT , OUTPUT , TAPE6=QUTPUT > 

190  PARAMETER (ND I M  =70) 

200  PARAMETER <ND I M2  =-NDIM) 

210  DIMENSION  CMOD ( ND I M ) , CARG ( ND IM > , VMOD ( NDI M) , VARG ( ND I M) 

220  COMPLEX  C ( ND I M ) , PS I ( ND IM2 : ND IM) , V ( ND I M ) , VA , COEFF , E , C0EFF2 

230  COMPLEX  R (ND IM ) , P ( ND IM ) , Q < ND I M) 

240  REAL  MU , L , K 

250  COMMON  P I , K , MU , EPS , W 

260  OPEN (UN I T=4 , F ILE=  CURRENT ' ) 

270  REWIND  4 

280  C 

290  C  DECLARE  CONSTANTS  EPS  AND  MU,  PERMITTIVITY  AND  PERMEABILITY 

300  C 

310  5  EPS  =  8.B4194E-12 

320  MU  -  1 . 256637E-06 

330  PI  -  3.14159 

340  C 

350  C  READ  IN  FREQUENCY,  LENGTH  AND  RADIUS  OF  THE  WIRE,  VOLTAGE 

360  C  THE  CENTER  OF  THE  WIRE,  AND  THE  NUMBER  OF  SEGMENTS. 

370  C 

380  10  WRITE (6,900) 

390  READ  *,  F 

400  WRITE (6,905) 

410  READ  *,  L . A 

420  WRI TE (6,910) 

430  15  READ  *,VA 

440  WR I TE (6,915) 

450  READ  *,  N 

460  WRITE (6,916) 

470  916  FORMA r ( 2X, ' INPUT  FAR-FIELD  DISTANCE',/) 

480  READ  * , RO 

490  C 

500  C  CALCULATE  ANGULAR  FREQUENCY,  WAVE  NUMBER,  SEGMENT  LENGTH 

510  C 

520  N1  -  N*(N+2> 

530  20  W  =  2*PI*F 

540  K-W*SQRT ( MU*EPS ) 

550  DL  =  L/ (FLOAT (N)  ♦  1.) 

560  C 

570  C  ASSIGN  VOLTAGE  VECTOR 

580  C 

590  CALL  VOLTS <N,V,VA) 

600  C 
610  C 

620  C  INITIALIZE  CURRENT  VECTOR 

630  DO  25  I = 1 , N 

640  C ( I )  =  (O. ,0. ) 

650  25  CONTINUE 

660  C  ASSIGN  IMPEDANCE  MATRIX 

670  C 

680  CALL  IN I TPS I (N,PSI ,A,DL> 

690  C 
700  C 


CALCULATE  CURRENT  VECTOR 


710 

c 

720 

CALL  CQNGRAD  (N,C,V,R,F',Q,  PS  I  > 

730 

DO  30  I  =■  l,N 

740 

CMOD ( 1 )  -  CABS (C ( I ) ) 

750 

Cl  -  A I  MAG ( C ( I ) ) 

760 

CR  *  REAL ( C ( I ) ) 

770 

CARG(I)  -  AT AN2 (Cl, CR ) *180 /PI 

780 

WRITE (6, 920)  I  , CMOD < I )  , CARG ( I ) 

790 

30 

CONTINUE 

800 

WRITE  c 4,925 >N 

810 

925 

FORMAT (13) 

820 

WRITE (4,930) 

830 

930 

FORMAT  (  C,  MODULUS  OF  ) 

840 

WRITE (4,935) ( CMOD ( 1 ) ,CMOD(I) ,I=1,N> 

850 

935 

F ORMAT  <  50 ( E 1 6 . 8 , / )  ) 

860 

WRITE (4,940) 

870 

940 

FORMAT  <  ' C ,  ARGUMENT  OF ' ) 

880 

WRITE (4,945)  <  CAR6 ( I )  , CARG (I)  , 1-1 ,N> 

890 

945 

FORMAT (50 (F 10. 3,/) > 

900 

COEFF  -  (0. , 1. >*SQRT (MU/EPS) *CEXP( (0 

910 

DO  56  J-0,18 

920 

THETA  »  10.  *J 

930 

THETAR  »  PI/180. *THETA 

940 

E  =  (0. ,0. ) 

950 

PHI  -  K*COS (THETAR) 

960 

C0EFF2  -  TAN ( THETAR) *COEFF*SI N ( 0. ! 

970 

DO  57  I  -  1  ,N 

980 

R 1  -  DL* ( FLOAT ( I ) -0 . 5*FLQAT ( N+ 

990 

E  -  E  +  C(I) *C0S(PHI*R1> 

1000 

57 

CONTINUE 

1010 

E  -  COEFF2*E 

1020 

EMOD  -  CABS (E) 

1030 

WRITE (6,956) THETA, EMOD 

1 040 

WR I TE ( 4 , 956 ) THETA , EMOD 

1 050 

956 

FORMAT (2X,  THETA  =  ,F5.1,2X,'E  » 

1060 

56 

CONTINUE 

FORMAT (2X,  THETA  -  ,F5.1,2X,'E  »  ',E10.4> 

CONTINUE 

FORMAT (2X INPUT  FREQUENCY  IN  HERTZ ',/> 

FORMAT <2X , ■ INPUT  LENGTH  AND  RADIUS  OF  THE  WIRE',/) 

FORMAT <2X,  INPUT  COMPLEX  VOLTAGE  APPLIED  TO  THE  CENTER  OF  THE 


l  WIRE  IN  THE  FORM  <VR,VI)',/> 

FORMAT <2X,  INPUT  THE  NUMBER  OF  SEGMENTS 
FORMAT (2X, 1 3 , 5X , E 1 2 . 4 , 5X , F7.  2) 

STOP 

END 


ODD  I NTEGER '  ,  /  > 


SUBROUTINE  VOLTS  ASSIGNS  VOLTAGE  VECTOR 


SUBROUTINE  VOLTS(N,V, VA) 
COMPLEX  V ( N ) , VA 
I  »  0 

IF ( 2* l + 1 . EQ . N ) GO  TO  lO 
V  ( 1  +  1  >  -  (O. ,0. > 

V (N-I )  »  <0. ,0. ) 


30 


1  =  1*1 
SO  TO  5 
V  (  I  ♦ 1 ) =VA 
RETURN 
END 


C0N6RAD  ITERATES  TO  THE  SOLUTION  VECTOR  C 


SUBROUT I NE  CONGRAD (N,C,V,R,P,Q,PSI> 

COMPLEX  C(N) , V (N) ,R(N) ,P(N) ,0(N> ,PSI ( -N i N ) 

CALL  SQUARE <N,V,V2) 

CALL  OPERAT (N,C,PSI ,0) 

J  -  0 

DO  100  1=1, N 

R ( I )  =  0(I>-V(I) 

CONTINUE 

CALL  SQUARE (N,R,R2) 

ERROR  =  R2/V2 

IF (ERROR. LT. IE-8) GO  TO  BOO 

CALL  COPERAT ( N , R , PS  I ,0) 

CALL  SQUARE <N,0,B> 

IF (B.LT. lE-OS)GO  TO  2 

B  =  1/B 

GO  TO  3 

B  =  l.E+OB 

CONTINUE 

DO  200  I  =  1  ,  N 

P<I>  =  — B*0 ( I ) 

CONTINUE 

CALL  OPERAT <N,P,PSI,0) 

CALL  SQUARE <N,0,T> 

IF  <  T. Lr . IE-8) GO  TO  3 

T  «  l/T 

GO  TO  6 

T  =  l.E+OB 

MR I TE <6, 900 >  J, ERROR 

FORMAT (2X,  '  ITERATION' , 13, 2X,  ERROR  =  '.E10.4  ) 
J  -  J*1 
DO  300  1=1 ,N 

C  <  I )  =  C ( I ) *T*P  < I ) 

CONTINUE 

CALL  OPERAT  <N,P,PSI ,0) 

DO  400  1=1 ,N 

R  ( I )  =  R ( I ) *T*0< I ) 

CONTINUE 

CALL  OPERAT (N,C,PSI ,0) 

DO  15  1=1, N 

0(1)  =0(1)  -  V(I) 

CONTINUE 

CALL  SQUARE < N , 0 , ERROR ) 

ERROR  =  ERROR /V2 

IF (ERROR. LT . IE-8) GO  TO  800 

CALL  COPERAT (N,R,PSI ,0) 


1810 

CALL  SQUARE (N, 0,8) 

1820 

IF  (B.LT. IE-8 ) GO  TO  18 

1 830 

8  -  1/8 

1840 

GO  TO  19 

1850 

18 

8  =  l.E+08 

1860 

19 

CONTINUE 

1870 

DO  500  I  - 1  ,  N 

1880 

P  (  I  )  *  P ( I > -8*0 (  I  ) 

1890 

500 

CONTINUE 

1 900 

CALL  OPERAT (N,P,PSI ,0) 

1910 

CALL  SQUARE (N,0,T> 

1920 

IF (T.LT. IE-8) GO  TO  16 

1930 

T  =  1/T 

1940 

GO  TO  17 

1950 

16 

T  =  l.E+O  1960  17  WRITE (6,900) J , ERROR 

1970 

7F<J.LT. 100)80  TO  1 

1980 

WRITE (6,930) 

1990 

930 

FORMAT (2X,  FAILED  TO  CONVERGE  AFTER  100  ITERATIONS 

2000 

GO  TO  801 

2010 

800 

CONT I NUE 

2020 

WR I TE  <  6 , 900 ) J , ERROR 

2030 

WRITE (6,920) 

2040 

920 

FORMAT <2X,  THE  ANSWER  IS) 

2050 

801 

RETURN 

2060 

END 

2070 

SUBROUT  I NE  GAUSS ( FCN , D , R , U , N , PS I ) 

2080  C 

2090  C  96  POINT  GAUSSIAN  RULE  INTEGRATION 

2100  C 

2110  COMPLEX  Z,PSI ,FCN,PSI 1 ,PSI2 

2120  DIMENSION! 48)  ,W(48) 

2130  DATA  A/ . 999689SOEOO ,  . 99836438E00 ,  . 99598 184E00 ,  . 99254390EOO 

2>4o  1  . 988054 13E  00,  .98251 726E  00,  . 975939 17E  00,  . 96832683E  00, 

2150  1  . 95968829E  00,  .95003272E  00,  .93937034E  OO,  .92771246E  00, 

2160  1  . 91507 142E  00,  .90146064E  00,  . 88689452E  00,  .87 13885 IE  00, 

2170  1  . 854959C3E  00,  .8376235 IE  00,  .81 94003 IE  00,  . 80030874E  00, 

2180  1  . 78036904E  00,  . 75960234E  00,  . 73803064E  00,  .71567681E  OO, 

2190  1  . 69256454E  00,  .6687 183 IE  OO,  .64416340E  OO,  . 61B925B4E  00, 

2200  1  . 59303236E  OO,  .56651042E  00,  .5393881  IE  00,  .51169418E  00, 

2210  1  . 48345797E  00,  .45470942E  00,  . 42547899E  00,  . 39579765E  00, 

2220  1  . 36569686E  00,  . 33520852E  OO,  . 30436494E  00,  .273 1988 IE  00, 

2230  1  . 24 1 743 1 6E  00,  .210O3131E  00,  .17B09688E  00,  .14597371E  00, 

2240  1  .11 369585E  00,  . 8 1 297495E-0 1 ,  . 488 1 2985E-0 1 ,  . 16276745E-01  / 

2250  DATA  W/.79679207E-03,. 18539608E-02, . 29 10731 8E-02, . 39645543E-0 

2260  1  . 50 1 42027E-02 ,  . 60585455E-02 ,  . 70964708E-02 ,  . 8 1 268769E-02 , 

2270  1  .91 48671 2E— 02 ,  . 10160771E-01 ,  . 1 1 162102E-01 ,  . 12151605E-01 , 

2280  1  . 13128230E-01 ,  . 1 4090942E-01 ,  . 1503872 IE-01 ,  . 15970563E-0 1 , 

2290  1  . 16885480E-01 ,  . 1 7782502E-01 ,  . 1 8660680E-01 ,  . 19519081E-01 , 

2300  1  . 20356797E-01 ,  . 2 11 72940E-0 1 ,  . 2 1 966644E-0 1 ,  . 22737070E-0 1 , 

2310  1  . 23483399E-01 ,  . 24204B42E-O1 ,  . 24900633E-01 .  . 25570036E-0 1 , 

2320  1  . 262 1 234 1E-0 1 ,  . 26826867E-01 ,  . 274 1 2963E-0 1 ,  . 27970008E-0 1 , 

2330  1  . 2849741 lE-Ol ,  . 289946 i 4E-01 ,  . 2946 1 090E-0 1 ,  . 29896344E-0 1 , 

2340  1  . 302999 15E-01 ,  . 306 7 1 376E-0 1 ,  . 3 1 0 1 0333E-0 1 ,  . 3 1 3 1 6426E-0 1 , 

2350  1  . 31589331E-01 ,  . 3 1 828759E-0 1 ,  . 32034456E-0 1 ,  . 32206205E-0 1 , 


32447164E-01 


725161 19E-01 


325506 14E-01/ 


2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

25 1 0 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2600 

2610 

2620 

2630 

2640 

2650  C 

2660  C 

2670  C 

2680  C 

2690  C 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

28 1 0 

2820 

2830 

2840  C 

2850  C 

2860  C 

2870  C 

2880  C 

2890 

2900 


1  •  -* — -  _'4  ..82  iE-O  If  • 

NPOINT  -  48 
PS  I  »  ( 0 .  ,  O .  ) 

G  -  U/2. /FLOAT  (N) 

DO  2  I  =  1  ,  N 
Z  =  <0.  ,0.  > 

DO  1  K  =  1, NPOINT 
ARGA  =  G*(A(k>  +  2*1-1.) 

ARGB  =  G* ( -A ( K )  +  2*1-1.) 

Z  =»  Z  +  (FCN  <  ARGA ,  D ,  R  >  +  FCN ( ARGB , D , R )  ) *W <  ►  ) 

1  CONTINUE 

PSI  -  PSI  +  Z 

2  CONTINUE 
PSI1  -  G*PSI 
PSI  =  <0.,0.) 

G  =  -G 

DO  4  I  »  1  ,  N 
Z  ■  <0. ,0. ) 

DO  3  K  *  l.NPOINT 
ARGA  -  G* (AIK)  +  2*1-1.) 

ARGB  »  G* ( —A ( K  >  +  2*1-1.) 

Z  -  Z  +  (FCN (ARGA . D , R )  +  FCN (ARGB , D , R >) *W <K ) 

3  CONTINUE 

PSI  -  PSI  ♦  Z 

4  CONTINUE 

PS I 2  -  G*PS I 

PSI  -  PSI 1  -  PS  1 2 

RETURN 

END 


INI  TPS I  CALCULATES  THE  PSI  VALUES 


SUBROUTINE  IN I TPS I ( N , PS I , A , DL ) 
COMPLEX  PSI <-N«N) ,CINT , GFUNCT ,F 
REAL  K.MU 

E  X  TERNAL  C I NT , AFUNCT , BFUNC T , GF  UNCT 
COMMON  PI , K.MU, EPS, U 
U*0. 5*DL 

F  *  DL/ ((0. ,-l.) *W*4. *P I *EPS> 

DO  1  I—  N.N 
D  -  I *DL 

CALL  GAUSS  (CINT.D,  A.U,  1  ,F'SI  ( I )  > 
PSI ( I )  »  F*PS I ( I ) 

1  CONTINUE 
RETURN 
END 


CINT  CALCULATES  THE  COMPLEX  INTEGRAND 


COMPLEX  FUNCTION  C I NT ( ARG , D 1 , A 1 > 
COMPLEX  GFUNCT, B.D 


2910 

REAL  k 

2920 

COMMON  P I  ,  K 

2930 

A  ”  k**2 

2940 

B  »  (O. ,-l. )*K/ (BFUNCT (ARG, D1 ,A1>**. 5) 

2950 

C  -  -(1  +  (K**2) *AFUNCT( ARG, D1 >) /BFUNCT (ARG, D1 

,A1 ) 

2960 

D  «  3*(0. , 1. >*K*AFUNCT (ARG, Dl) 7 (BFUNCT (ARG, D1 , 

A1 ) ** 1 

2970 

E  =■  3*AFUNCT  (ARG, D1 )/  (BFUNCT  (ARG, D1  ,A1 )  **2> 

2980 

C I NT  =  *A  +  B  +  C  +  D  +  E >  *GFUNCT ( ARG , D 1 , A 1 ) 

2990 

RETURN 

3000 

END 

3010  C 
3020  C 

3030  FUNCTION  AFUNCT (ARG , D > 

3040  AFUNCT  -  CARG-D)**2 

3050  RETURN 

3060  END 

3070  C 

3090  C 

3090  FUNCTION  BFUNCT (ARG , D , A 1 ) 

3100  BFUNCT  »  AFUNCT (ARG, D)  +  Al**2 

3110  RETURN 

3120  END 

3130  C 

3140  C 

3150  C  GFUNCT  CALCULATES  THE  GREEMI/UNCT I  ON 

3160  C 
3170  C 

3180  COMPLEX  FUNCTION  GFUNCT (ARG, D ,A1 > 

3190  REAL  K 

3200  COMMON  P I , K 

3210  R  »  SORT ( BFUNCT ( ARG, D,A1) > 

3220  GFUNCT  -  CEXP( (0. ,-l. >*K*R)/R 

3230  RETURN 

3240  END 

3250  C 

3260  C 

3270  C  COPERAT  IS  THE  ADJOINT  OPERATOR 

3280  C 
3290  C 

3300  SUBROUTINE  COPERAT (N, R , PS I ,0> 

3310  COMPLEX  PS I ( -N ;N) ,R(N) ,Q(N> 

3320  1  CONTINUE 

3330  DO  200  I  *  1 , N 

3340  0(1)  «  (O. ,0. ) 

3350  DO  100  J«1,N 

3360  M  -I-J 

3370  0(1)  »  0(1)  +  CONJG (PS I (M) >  *R ( J  > 

3380  1 00  CONT I NUE 

3390  200  CONTINUE 

3400  RETURN 

3410  END 

3420  C 
3430  C 
3440  C 
3450  C 


SQUARE  MULTIPLIES  A  VECTOR  BV  ITS  COMPLEX  CONJUGATE 


LISTING  NO.  2 


1  00  C 

110  L  PROGRAM  WhN  I  r D  CALCULATES  THE  LURRENTS  UN  A  WIRE  aNTENNA 

120  C 

120  C  WANTED  CALCULATES  ! HE  CURRENTS  ON  A  WIRE  ANTENNA  WITH  THE 

140  C  SOURCE  AT  THE  CENTER.  ‘C>  13  THE  CURRENT  VECTOR,  'V)  IS  THE 

150  C  APPLIED  VOLTAGE  VECTOR  AND  OF'ERAT  IS  THE  INTEGRAL 

160  C  OPERATOR  SUBROUTINE.  WANTED  SOLVES  FOR  THE 

170  C  CURRENTS  USING  THE  CONJUGATE  GRADIENT  METHOD. 

130  PROGRAM  WANTED' I NPUT  ,  T«r  E5=  I NPUT  ,  OUTPUT  ,  TAPE6=GU  TF'UT  .■ 

190  PARAMETER l ND I M  =  ~Bh 

200  PARAMETER <ND I M2  =-NDIM) 

213  D I  MENS  I  ON  CMOD  (NDIM)  ,  C  ARb  t  ND  I M  /  ,  VMOD  <  ND  1 M .  .■  ■  ND  I M > 

220  COMP  LEX  C  i  ND  I M  >  ,  PS  I  ( ND  IM2:  ND  IM )  ,  V  -  ND 4  . -•  ■  .  '  ,  ,.OEF'F2 

270  COMPLEX  R  (NDIMf  ,p  •  NDIM)  ,0  . 2 ■ NDIM2: ND I M ■ 


240 

REAL  MU,L,P 

250 

COMMON  PI,. 

■  •  ,  w 

260 

OPEN  •  ■  ■  ■  ’  -  . 

iL  CUR  PEN  T 

270 

T  ■_  )  ,  •  •> 

230  C 

290  •.  DLCLaRE  CONSTANTS  LRS  AND  MU,  PERMIT  I  IVI  I 

300  C 

310  5  EPS  -  0.fJ4f'4L  12 

320  MU  -  1 . 25oo37E- 06 

770  T I  =  3.14159 

740  C 

750  -  RLmD  IN  ERLUUENUX  .  LENGTH  uND  RaDIUs  U) 

7j0  c  the  center  or  ’•ml  wire,  «no  the  number 

”0  L 


:co 

I  V.' 

w>- :  * c 

’-J  ,  vuvJ 

.'•’0 

'Lt-I.- 

* 

,  r 

>00 

WR  ;  1  L 

•>,  .•O''. 

1  :  X 

•-•LhD 

- 

,  L  ,H 

■>:> 

WF  1  It 

•  - 

■-  ,  -7  1  U 

4  0 

L 

'  'Li  *L' 

«• 

,  Vh 

>40 

i  L 

.j  ,  O  l  \ 

•ISO 

*  C  J-4  0 

♦ 

,  n 

>  0 

i  ri. 

j ,  •  i  j 

0  '  .  if»  •  M  ILl.D  I'l'PET.l.  ,  • 

;  o  m-;ad  *,>t 

Lri  '  T  ftN  ,'jL  f-t:  :  '  L.  U '  j f.  V ).  Jljl't... 

N  J 

w  c  »t-  ;  *  ■ 

•  -  W  <  OR  1  .  M1  '  • 

CL  L  -rl  j.,'  i;  ;  . 

•  "  .  '■  aGL  VECTOR 

ALT.  .-’UL  T S  ■  N  ,  V  ,  .  A ) 


INI’ IALI2E  CURRENT  VECTcR 
DO  25  I  *  1  .  N 

L  1  i  ’  :  )  0  .  ,  0  . 

CONTINUE 

ASSIGN  IMPEDANCE  MaTRI> 

CALL  IN I  TPS  I  IN. PC  I , m , ul ) 

CALCULATE  CURRENT  VECTOR 


.  J 

1 3  ..  U 

■-•it) 


■  .’0  ■_ 
j  v 1 
..■1)0 
'  ' '  c 

62  0  L 
6  70 
6  40 
'.50 
6o0  T. 
6  ’0  C 
t.30 
6^>0  C 
T00  C 


t  AND  RERMEhBILI  i  i 


THE  WIRE,  VOL r AGE 
OF  SEGMENTS. 


.-GHENT  LENGTH 


37 


-10 

C 

770 

lALL  LONG! (AD  ■  N ,  C  ,  V  ,  R  ,  P  ,  U ,  Z  • 

750 

DO  30  I  =  I  .  N 

-40 

C MOD  <  I  ‘  =  CABS  i.  L  ‘  I  ■>  / 

750 

Cl  =  A I MAG  t C '  I  '  ' 

760 

CR  =  REAL  '  C  (  I  >  ) 

7  70 

CARG < I >  =  A  T AN7 'Cl , CR ) *  1 80 / P I 

780 

WR  I  TE  1 6 , 920 I  ,  CMOD  ■  1  ,  CARG  <  I  .■ 

790 

50 

CONTINUE 

900 

WR I TE  <  4 , 975 ) N 

310 

975 

FORMAT (15) 

870 

WRITE (4, 930) 

350 

950 

FORMAT ( • C .  MODULUS  OF  > 

340 

WRITE (4,935)  t LMQD  < I )  ,CMOD( I )  , 1  =  1  ,N 

850  935  FORMAT i50 (Elib. 3, /> > 

3o0  WRITE (4,940) 

370  940  FORMAT «C,  ARGUMENT  OF  : 


980  WR I TE  <  4 , 945 >  •  CARG  <  1 ,  CARG  <  I )  ,  I  *  1  ,  N  > 

890  945  FORMAT  i50 '  F10.  5,  /  • 

900  COEFF  =  (0.  ,  t .  •  *SQRT  iMU/EPS)  *CExP  i  <0.  .-1 .  >  **  *RO>  >  *PI*R0) 

910  DO  5d  J=0,18 

970  THETA  =  10.  *J 

950  THETAR  =  PI/ 130. ♦THETA 

940  E  -  i0. ,0. ) 

950  PHI  =  \  *C03 ( THETAR) 

960  COEFF7  -  TAN  <  THETAR' )  *CQEFF*SIN  ( 0. 5*DL*PH I > 

970  DO  57  I  =  1 ,N 

OQ0  R1  =  DL*  (FLOAT  <  I  >  -0.  5*FLQAT  iN+  1  ! 

990  E  =  E  +  C < I ) *CQS (PHI*R1 ) 

1000  57  CONTINUE 

1010  E  =  COEFF2*E 

1070  EMOD  =  CABS ( C ) 

1050  WRITE <6. 959) THETA, EMOD 

1040  WRITE (4,959) THETA, EMOD 

1050  956  FORMAT! 7a,  THETA  =  ,F5.l,7X,  E  =  ,E10.4) 

1060  56  CONTINUE 

1070  900  FORMAT ( 7X  ,  INFUT  FREQUENCY  IN  HERTZ  ,/> 

1080  905  FORMAT ' 7X,  INPUT  LENGTH  AND  RADIUS  OF  THE  WIRE  ,/> 

1090  910  FORMAT • 2X ,  INFUT  COMPLEX  VOLTAGE  APPLIED  TO  THE  CENTER  OF  THE', 

1100  1  WIRE  IN  THE  FORM  <VF*,Vl.'  ,  /  > 

1110  915  FORMAT 17X,  INPUT  THE  NUMBER  OF  SEGMENTS  =  ODD  INTEGER  ,/i 
1170  9 70  FORMAT ‘ 7a , I5,5x ,E17. 4,5a ,F7. 7) 

1150  STOP 

1  1  40  END 

1  150  C 
l  1  C 

1170  C  SUBROUTINE  VOLTS  ASSIGNS  VOLTAGE  VECTOR 

1180  C 
1  190  C 

1700  SUBROUTINE  VOLTS N , 'V ,  VA ) 

1710  COMFLEA  V  ( N ;  .  VA 

1770  I  =  0 

1750  5  IF  (7*1*1 .  EQ.N.'GO  TO  10 

1  740  V  (  I  >  1  •  =  i  0 .  ,  0 . 

1750  V (N-I  >  =  <0.  ,  0.  ' 


V  -V  V  7  > 


mxuL. 


1260 

1  =  1*1 

1  270 

GO  TO  5 

1  230 

10 

V  (  I  ♦  1  >  =  v'A 

1  290 

RETURN 

1  200 

END 

1210 

f* 

1  320 

c 

1  330 

c 

CQN6RAD  I T  Ef-ViTES  TO  THE  SOLUTION  VECTOR  L 

1340 

c 

1350 

c 

1360 

SUBRGUT  I NE  CONOR  AD  •,  N ,  C  ,  V  .  R  ,  R  ,  U ,  2  ) 

1370 

COMF'LE  C  IN)  ,  V  1  N  '  ,PiN(  ,  P  i  N  ,  u  ( N )  ,  Z  ■'  -N :  N ) 

1380 

CALL  SQUARE . N , V , V2 • 

i  3^0 

CALL  OPERA F ■ N.C.Z.O' 

1  400 

J  =  0 

1410 

DO  100  1=1 , N 

1420 

R  I  1  =  /  • 

1430 

100 

CONTINUE 

1440 

CALL  SQUARE !N,R, R 2: 

1450 

ERROR  =  R2/V2 

1460 

IF (ERROR.LT. IE -9100  TO  800 

1470 

CALL  COPERAT  « N ,R , Z , 0) 

1480 

CALL  SQUARE (N,G,B> 

1490 

IF ( B. LT. IE-00; GO  TO  2 

1500 

B  *  l/B 

1510 

GO  TO  3 

1520 

—1 

B  *  1.L-C08 

1530 

CONT INUE 

1540 

DO  200  I=1,N 

1550 

PFI.'  =  -B  *0  <  I  ■ 

1560 

200 

CONTINUE 

1570 

CALL  OF'ERAT  <  N ,  P  ,  Z  ,  0 ) 

1580 

CALL  SQUARE  IN,  0,  T> 

15Q0 

IF  \  T  LT.  1 E -8 GO  TO  5 

.  600 

T  =  1/T 

1610 

GO  TO  6 

1620 

5 

T  =  I.E+08 

1630 

l 

WP I TE  <  6 , 900 ) J , ERROR 

1640 

900 

FORMAT (2X  ,  ITERATION  , I’,2X,  ERROR  =  ,E10.4 

.650 

1 

J  =  J*1 

1660 

DO  300  I  =  1  ,  N 

1670 

c  a  :•  =*  c  <  i  j  «-t*p  ■  i  .• 

lo80 

300 

CONTINUE 

.690 

caul  opera  r ( n . p . z , o  > 

1700 

DO  400  I = 1 „N 

’10 

RT)  =  R  a  >  +-T  *0  (  I  > 

720 

400 

CONTINUE 

7  30 

CALL  OPEPAT  (N.'_  ,  Z  .0  ■ 

340 

DO  15  I = 1 , N 

1750 

0  <  1 >  =0(1.'  -  V  <  i  > 

’60 

15 

CONTINUE 

770 

CALL  SQUARE 'N.O. ERROR.' 

L  780 

ERROR  =  ERROR.-  V 2 

790 

IFvERROP.LT. IE-8)00  TO  300 

800 

CALL  COPERAT (N,R, Z, O' 

810 

320 

830 

840 

350 

18 

8o0 

1  9 

8  70 

380 

390 

500 

900 

910 

920 

930 

940 

950 

16 

9  60 

17 

970 

980 

9  90 

930 

000 

010 

300 

020 

030 

040 

920 

050 

801 

060 

070 

080 

C 

090 

C 

100 

C 

1  10 

120 

1  30 

140 

150 

lo0 

170 

190 

190 

200 

210 

220 

230 

240 

250 

CALL  SQUARE  •  f J .  0  .  B  • 
ir  iu.lt.  1C.-B.  go  ro  ic 
b  =  i .  b 

CO  TO  19 

&  =  i.e+oq 

CONTINUE 
DO  500  I  =  1  ,  N 

F  •  I  >  —  K  ■  i  '  b  *  1 J  I 

CONTINUE 

CALL  OPERA T N,F , £ ,u • 

CALL  SQUARE  •  N  ,  0  .  '  .■ 

I F  (  T  .  L  T  .  1  E  -  S  )  ijO  TO  io 
T  =  1/T 
00  TO  17 
T  =  l.E+08 
WRITE<o,«00> J .ERROR 
IFiJ.LT.  100,'GO  TO  1 
WRITE  i.«j,930.' 

FORMAT <2X ,  FAILED  TO  CONVERGE  AFTER  100  ITERATIONS  > 

GO  TO  301 
CONTINUE 

WR I TE  <  6 , °O0 ) J . ERROR 
WR I TE  (6 , 920 

FORMAT (OX,  THE  ANSWER  IS  > 

RETURN 

END 

SUBROUTINE  GAUSS <FCN,D,R,U,N, RSI) 

96  POINT  GAUSSIAN  RULE  INTEGRATION 

COMPLEX  Z , PS  I , FCN , PS  I  1 , PS 1 2 
DIMENSION  A ( 48 )  ,  W ' 48 ) 

DATA  A/.  99968950E00,  . 99836438E00 ,  .  99593 1 84E00 ,  .  992543°0E00 


1 

. 988054 13E 

00, 

.98251726E 

00, 

. 975939 17E 

00, 

. 96832683E 

00, 

1 

.95968829E 

00, 

. 95003272E 

00, 

. 939  37034E 

00, 

. 9277 1 246E 

00, 

1 

.91507142E 

00, 

. 90 1 46064E 

00, 

.B8689452E 

00, 

. B7 1 3885 IE 

00, 

1 

. 85495903E 

00, 

. 83762351E 

00, 

.81 94003 IE 

00, 

.80030874E 

00, 

1 

. 78036904E 

00, 

. 75V60234E 

00. 

. 73803064E 

00, 

. 71567681E 

00, 

1 

. o9256454E 

00, 

. 0687 183 IE 

00, 

. 64416340E 

00, 

.61 892584E 

00, 

1 

. 59303236E 

00, 

.56651042E 

00, 

.5393881  IE 

00, 

. 51 169418E 

00, 

1 

. 48345797E 

00, 

. 45470942E 

00, 

. 42547899E 

00, 

. 39579765E 

00, 

1 

. 36569636E 

00, 

. 33520852E 

00, 

. 30436494E 

00. 

. 27319881E 

00, 

1 

. 2417 43 1 6E 

00, 

. 2 1 003 1 3 1 E 

00, 

.  1  ‘7S09688E 

00, 

.  1 459737  IE 

00, 

l 

. 1 1369585E 

00, 

.81 297  495E  • 

-01, 

. 488 1 2985E- 

-01  , 

.  16276  745E- 

-01  / 

DATA  W/ . 79679207E-03.. 1 853960SE-02 , . 29 1 073 1BE-02 , . 39645543E-0 


260 

1 

. 50142027E -02 , 

. O0585455E-02, 

. 70964708E-02, 

.81 2687 cj9E-02  , 

270 

1 

. 9148671 2E-02, 

. 10160771E -01 , 

. 1 1 162102E-01 , 

. 12151 605E-01 , 

280 

1 

. 1 3128230E-01 , 

.  1 4090942E  -01 , 

. 15038721E -01 , 

. l59705o3E-01 , 

290 

1 

. 16885480E-01 , 

. 17782502E-01 , 

. 18660680E-01 , 

.  1 95 1 908 1 E  -01 , 

300 

1 

. 20356777E-01 . 

. 21 172940E-01 , 

. 2196oo44E-01 , 

. 22737070E-01 . 

310 

i 

. 23483399E-01 , 

. 24204842E-01 , 

. 249006 33E-01 , 

. 25570036E-01 , 

320 

1 

. 26212341E-01 , 

.2682O867E-01 , 

. 27412963E-01 , 

. 27970008E-01 , 

330 

1 

. 2849741  IE-01 , 

. 289946 14E-01 . 

. 2°4o 1 090E-0 1 , 

. 29896 344E-01 , 

340 

1 

. 302999 15E  01 , 

. 3067 1 37 6E-0 1  , 

. 31010333E-01 , 

. 31316426E-01 , 

350 

1 

. 3 158933 IE-01 , 

. 31828 759E-01 , 

. 320  3445oE-0 1 , 

. 32206205E-01 , 

40 


1 


-•1:447  X 64L-  01 


.'■-jI'jI'.vL  01 


325506 14E-01/ 


1360 
1370 
1380 
1390 
1400 
1410 
1410 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1510 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1610 
1630 
1640 
1650  C 
1660  C 
1670  C 
1680  C 
1690  C 
1700 
1710 
1~10 
1730 


. 31343813L-01 ,  . 

NFOINT  -  48 
PSI  =  '0. ,0. ) 
a  =  u/i.  ,floht  >n.' 

DU  1  l  =  1  ,  N 
Z  =  10. ,0. > 

DO  L  1  =  l.NPUINT 

AFGA  =  G  *  '  A  >  1  t-  1  *  £  - 1 .  ) 

ARGB  =  ••  +  !■»  I  -  l .  ,i 

Z  =  Z  *  <FCN>hRgA,D,R)  *  FCN  ■ ARGB .  D ,  R  >  j  *W  <  1  > 

CONTINUE 

PS  I  =  PS  I  Z 

CONTINUE 

PS  II  =  G*F'S  I 

PSI  ■=  .0. ,0. ) 

a  =  -a 

DO  4  I  -  1  ,  N 
Z  =  <  0 .  ,  0 .  ) 

DO  3  k  =  l,NPOINT 

AP'GA  =  G*  ( A  ( k )  ♦  1*1-1.) 

ARGB  =  G* ( -  A ( k )  +  1* I  —  1 .  > 

Z  -  Z  *  (FCN  ( ARGA ,  D  .  R )  +  FCN  ( ARGB ,  D ,  R )  )  *W  K ) 

CONTINUE 

PSI  =  PSI  *  Z 

CONT INUE 

PS  II  =  G*FSI 

PSI  =  PS  II  PS  II 

RETURN 

END 


INI  TPS  I  CALCULATES  THE  PSI  VALUES 


SUBRQUT INE  INITFSI (N , PS I , A , DL , Z > 

COMPLEX  PSI l  -N:N> , C INT , PS I A , PG I B , PS IC , PS I D , PS  IE . Z < -N: N) 
REAL  1 , MU 
EXTERNAL  CINT 


COMMON  F I ,k,MU,EPS,W 
U-0. 5*DL 


1--0 

DO 

1  I=-N,N 

1780 

D  -  I *DL 

!T9Q 

CALL  GAUSS  <  C  I  NT  ,  D ,  A ,  U ,  1  ,  P'S  1  <  I 

1810 

1  CONTINUE 

181  1 

DO 

100  J  =  1 

,  N 

1813 

DO  100  I 

=  1  ,N 

1314 

M  =  1- 

J 

1315 

P3IA  = 

PS  I  < M  >  • DL 

18 1  o 

PS  I B  - 

PSlH 

1817 

M  =  I 

J+1 

1813 

PS  I C  = 

PSI  iM)  .  DL 

1319 

M  =  I  - 

J-l 

1010 

PS  ID  = 

PS  I  ( M  >  DL 

1011 

PS  IE  - 

PSI  A 

1811 

M  = 

I  -  J 

41 


cac-7  C'M>  -  ■  a.  ,  1  .•  •  w«mu*-L'L  ♦re *•-. 

2829  Z<»>  -  2  •  M  •  *  1  .  W.-  •  Id.  .  1  .  •  U.F--'i •  •  F  b  I B  F  iiC  -FSXD^PGIE* 

282°  100  CUM  r  I  NULL 

2820  200  CONTINUE 

2831  f-'ET'JFN 

2832  END 
2840  C 

2850  L 

28o0  L  CIM!  CALCULATES  'I  HE  LUMPLE'  INrLGRwNU 

28  70  C 
2880  C 

2890  COMPLEX  FUNCTION  C I NT i ARG . D , o : 


23°  1 

PEAL  F 

28^2 

COMMON  F'l.F 

2893 

P  -  SOP:  T  '  A'*'* 

2  ♦  •. ARQ-D>  ♦♦2.' 

20°4 

LINT  =  CEXFhi 

0 .  ,  -  1  .  *F  ♦R  >  .•  »  4  *  K  I  •*  R  .• 

2095 

RETURN 

28^3 

END 

3250 

c 

3260 

c 

3270 

c 

CORERaT 

12  THL  hDJOINT  OF  th'HTQf' 

3280 

r 

3300 

SUBROUTINE  CORERAT •N,R,2,0' 

3310 

COMPLEX  P«N) 

,Q(N)  . Z • -Ns  N ' 

3  330 

DO  200  1  =  1 

,N 

3  340 

0  •  I  )  =  <  0 

.  .0.  > 

3  350 

DO  100  J 

-  1  ,N 

3  3o0 

M  =  I- 

J 

•-  -0 

0  •'  I  ■  = 

0(1'  +  lUNJG i i iM;  , *P  < j  • 

3330 

100 

CONTINUE 

3  390 

200 

CONTINUE 

3400 

RETURN 

3410 

END 

3420  C 
3430 

3440  C  SQUARE  MULTIPLIES  a  VEC 1  UP  Br  lib  COMPLEX  CONJUGATE 

3  i5U  C 
~  480  C 

34~0  SUBROUF1NC  SQUARE  i  N ,  G 1  , U2 ; 

3480  COMF LE  <  01  ■>!.' 

34 ?0  02=0. 

3500  C'O  100  I  =  1,N 

3510  02  =  02  ♦  Ql  i  I  >  *CUNJG  Ujl  (  1  >  .« 

3520  100  CONTINUE 

7530  PETUF<N 

3540  END 

3550  C 
35o0  C 

"3570  C  OPERA  T  is  THE  INTEGRAL  OPERA  TUB 

3580  C 
3590  C 

-.300  SUBRO'JT  I  NE  UF'ER  hi  ■  N  ,  I,  .  2 , 0  ’ 

33 1 0  COMPLEX  C  •  N»  ,0 • N;  , Z <  -N: N- 

3620  DO  200  1=1 ,N 

3330  0 ' I ■  =  • 0. ,0. ) 


DO  1 00  J  =  1 . N 
M  =  I  J 

o  .■  n  -  u  ■  i )  *  -  m  ■  *  o  ■  J  ’ 

100  CONTINUE 

200  CONTINUE 
RETURN 
END 


3o40 
3o50 
3  300 
36  70 
3680 
3o90 
3700 


42 


LISTING  NO.  3 


100  c 

110  C  PROGRAM  WANTMOM  CALCULATES  THE  CURRENTS  ON  A  WIRE  ANTENNA 

120  C 

130  C  WANTMOM  CALCULATES  THE  CURRENTS  ON  A  WIRE  ANTENNA  WITH  THE 

140  C  SOURCE  AT  THE  CENTER.  CZJ  IS  THE  IMPEDANCE  MATRIX,  (C)  IS 

150  C  THE  CURRENT  VECTOR,  AND  (V)  IS  THE  APPLIED  VOLTAGE  VECTOR. 

160  C  WANTMOM  SOLVES  FOR  THE  CURRENTS  USING  THE  METHOD  OF  MOMENTS. 

170  C 

180  PROGRAM  WANTMOM ( INPUT , TAPES- I NPUT , OUTPUT , TAPE6-0UTPUT > 

190  PARAMETER (NDIM  -  70) 

200  PARAMETER (ND I M2  -  ND IM**2 , ND IM3  -  -NDIM.NDIM4  -  NDIM* (NDIM+2) > 

210  D  f MENS  I ON  CMOD (ND IM)  , CARG <ND IM)  , VMOD ( NDIM)  , VARG <NDIM) 

220  COMPLEX  Z (NDIM2)  , C (NDIM)  , PS  I (NDIM3: NDIM)  , V (NDIM)  , VA 

230  COMPLEX  WA (ND I M4) , WK (ND IM > , COEFF , C0EFF2 , E 

240  REAL  MU , L , K 

250  COMMON  P I , K , MU , EPS , W 

260  0PEN(UNIT=4, FILE-  CURRENT ' > 

270  REWIND  4 

280  C 

290  C  DECLARE  CONSTANTS  EPS  AND  MU,  PERMITTIVITY  AND  PERMEABILITY 

300  C 

310  5  EPS  -  8. 84194E-12 

320  MU  -  1.256637E-06 

330  PI  -  3.14159 

340  C 

350  C  READ  IN  FREQUENCY,  LENGTH  AND  RADIUS  OF  THE  WIRE,  VOLTAGE 

360  C  THE  CENTER  OF  THE  WIRE,  AND  THE  NUMBER  OF  SEGMENTS. 

370  C 

380  10  WRITE  (6-, 900) 

390  READ  *,  F 

400  WR I TE ( 6 , 905 ) 

410  READ  *,  L, A 

420  WRITE (6,910) 

430  15  READ  *,VA 

440  WRITE (6,915) 

450  READ  *,  N 

460  ) 

470  916  FORMAT (2X ,' INPUT  FAR-FIELD  DISTANCE ' ,/> 

480  READ  *,RO 

490  C 

500  C  CALCULATE  ANGULAR  FREQUENCY,  WAVE  NUMBER,  SEGMENT  LENGTH 

510  C 

520  N1  -  N* (N+2) 

530  20  W  -  2*PI*F 

540  K«W*SQRT(MU*EPS> 

550  DL  -  L/ (FLOAT <N)  +  1.) 

560  C 

570  C  ASSIGN  VOLTAGE  VECTOR 

580  C 

590  CALL  VOLTS (N,V,VA> 

600  C 

610  C  ASSIGN  IMPEDANCE  MATRIX 

620  C 

630  CALL  I MP ( N ,  Z  ,  A , DL , PS I > 

640  C 

650  C  CALCULATE  CURRENT  VEC . OR 

660  C 

670  CALL  LEQS  (  Z  ,  N  ,  N 1  ,  V  ,  WA ,  WK ) 

600  DO  30  I  -  1,N 

690  CMOD ( I )  -  CABS ( V ( I ) ) 

700  Cl  -  A I MAG ( V ( I ) ) 


<iJT r.f 


7  10 

CR  =  REAL (V  C I )  ) 

720 

CARG(I)  =  ATAN2(CI ,CR) *180/PI 

730 

WRITE (6,920)  I , CMOD ( I ) , CARS < I ) 

740 

30 

CONTINUE 

750 

DO  40  I  ■=  1  ,  N 

760 

C  (  I )  =*  V  (  I ) 

770 

40 

CONTINUE 

780 

WRITE (4,925) N 

790 

925 

FORMAT (13) 

800 

WRITE (4,930) 

810 

930 

FORMAT (C,  MODULUS  OF  '  > 

820 

WR I TE ( 4 , 935 ) ( CMOD ( I ) , CMOD ( I ) ,1=1,1 

830 

935 

FORMAT (50 (E 16. 8,7) ) 

840 

WRITE ( 4 , 940 ) 

850 

940 

FORMAT!  C,  ARGUMENT  OF  ) 

860 

WR I TE ( 4 , 945 ) ( CARG ( I ) , CARG ( I ) ,1=1,1 

870 

945 

FORMAT (50 (F10. 3, / ) ) 

880 

COEFF  =  (O. , 1. >*SQRT (MU/EPS) *CEXP 

890 

DO  56  J=0, 18 

900 

THETA  =  10. *J 

910 

THETAR  =  PI/180. *THETA 

920 

E  -  (0. ,0. ) 

930 

PHI  »  K#COS( THETAR) 

940 

C0EFF2  =  TAN ( THETAR ) *COEFF*S IN 

950 

DO  57  I  =*  1  ,N 

960 

R1  -  DL* (FLOAT ( I ) -0. 5*FL0AT 

970 

E  »  E  +  C ( I ) *CQS ( PH I *R  l ) 

980 

5  7 

CONTINUE 

990 

E  =  C0EFF2*E 

1 OOO 

EMOD  =*  CABS  (E) 

1010 

WR I TE ( 6 , 956 ) THETA , EMOD 

20 

WR I TE ( 4 , 956 ) THETA , EMOD 

1030 

956 

FORMAT(2X,  THETA  -  ,F5.1,2X,'I 

1 040 

56 

CONTINUE 

1050  900 
1060  905 


FORMAT (2X,  THETA  -  ,F5.1,2X,'E  -  ,E10.4> 

CONTINUE 

FORMAT (2X, ' INPUT  FREQUENCY  IN  HERTZ',/) 

FORMAT <2X,  INPUT  LENGTH  AND  RADIUS  OF  THE  WIRE  ,/> 

FORMAT (2X ,' INPUT  COMPLEX  VOLTAGE  APPLIED  TO  THE  CENTER  OF  THE', 


1  WIRE  IN  THE  FORM  <VR,VI>  ,/) 

FORMAT <2X,  INPUT  THE  NUMBER  OF  SEGMENTS 
FORMAT (2X , I3,5X,E12.4,5X,F7.2) 

STOP 

END 


SUBROUTINE  VOLTS  ASSIGNS  VOLTAGE  VECTOR 


ODD  INTEGER  , /> 


SUBROUTINE  VOLTS (N , V , VA ) 
COMPLEX  V ( N ) , VA 
I  -  O 

IF2*I*1.E0.N)G0  TO  lo 

V  (  I  +  1  )  =■  ( 0 .  ,  0 .  ) 

V  < N- I )  »  (O.  ,0.  ) 

I  -  I  *  1 

GO  TO  5 
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1260  10  0(I+1)«0A 

1270  RETURN 

1 280  END 

1290  C 
1300  C 

1310  C  SUBROUTINE  IMP  ASSIGNS  IMPEDANCE  MATRIX 

1320  C 
1330  C 

1340  SUBROUTINE  I MP ( N , 2 , A , DL , PS 1 ) 

1350  COMPLEX  Z (N,N>  ,PSI <-NsN>  ,PS I A , PS  I B , PS IC , PS  I D , PSIE ,C INT 

1360  REAL  k , MU 

1370  EXTERNAL  Cl NT 

1380  COMMON  P I , K , MU , EPS , W 

1390  C 

1400  U  -  0. 5*DL 

1410  C  INTEGRATE  FROM  -U  TO  U  TO  GET  PSIS 

1420  DO  1  I  «  -N,N 

1430  D  =  I *DL 

1440  C 

1450  C  D  IS  THE  DISTANCE  FROM  INTEGRATING  POINT,  J  TO  OBVERSATION 

1460  C  POINT,  I. 

1470  C 

1480  CALL  GAUSS  (CINT,D,A,U,  1  ,  F'S  I  (I>  > 

1490  1  CONTINUE 

1500  DO  200  J«1,N 

1510  5  DO  100  I-1,N 

1520  M  -  I-J 

1530  PSIA  -  PSI (M) /DL 

1540  PSIB  -  PSIA 

1550  10  M  -  I-J+l 

1560  PSIC  =»  PSI  (M)  /  DL 

1570  M  -  I-J-l 

1580  PSID  »  PSI (M) /DL 

1590  PSIE  »  PSIA 

1600  15  Z(I,J>  -  <0.  ,  1 .  )  *W*MU*DL**2.  *F'SIA 

1610  Z(I,J)  -  Z(I,J)  +  1. /W/ <0. , 1 . > /EPS*(PSIB-PSIC-PSID+PSIE) 

1620  C  WRITE (6, 900)  I,J,Z(I,J> 

1630  100  CONTINUE 

1 640  200  CONT I NUE 

1650  900  FORMAT <2X, ' Z <  , 13, 13, ' )  =  .2E15.8) 

1660  RETURN 

1670  END 

1680  C 
1690  C 

1700  C  Cl NT  CALCULATES  THE  COMPLEX  INTEGRAND 

1710  C 
1720  C 

1730  COMPLEX  FUNCTION  C I NT l D , A , ARG ) 

1  740  REAL  K 

1  750  COMMON  PI ,k 

1760  R  -  SORT (A**2  ♦  <ARG-D)**2) 

1770  CINT  -  CEXP  <  (0.  ,-l .  )  *»■  *R)  /  (4*F  I*R) 

1780  RETURN 

1 790  END 

1800  C 


LINEAR  EQUATION  SOLUTION 


1810  C 
1820  C 

1830  SUBROUTINE  LEQS<Z ,N,N1 ,V,WA,WK> 

1840  COMPLEX  Z<N,N) ,V<N,1) , W A ( N 1 ) , WK ( N ) 

18S0  CALL  LEQ2C ( Z , N , N , V , 1 , N , O , WA , WK , I ER  > 

1860  RETURN 

1870  END 

1880  SUBROUTINE  GAUSS <FCN,D,R,U,N,PSI > 

1890  C 

1900  C  96  POINT  GAUSSIAN  RULE  INTEGRATION 

1910  C 

1920  COMPLEX  Z , PSI ,FCN, PSI 1 , PSI2 

1930  DIMENSION  A(48) ,W(48> 

1940  DATA  A/ . 9996S950E00 ,  . 99B3643BE00 ,  . 99398 1 84EOO ,  . 99254390E00 , 

1950  1  .988054 13E  00,  .98251726E  00,  . 975939 17E  00,  .96832683E  OO, 

1960  1  . 95968829E  OO,  .950O3272E  00,  .93937034E  00,  .92771246E  OO, 

1970  1  91507142E  00,  .90146064E  00,  . 886894S2E  00,  .871388S1E  OO, 

1980  1  . 85495903E  OO,  .8376235 IE  00,  . 81940031 E  00,  .80030874E  OO, 

1990  1  . 78O36904E  00,  . 75960234E  00,  . 73803064E  00,  .71S67681E  OO, 

2000  1  . 692S64S4E  00,  .6687 183 IE  00,  .64416340E  00,  .61892S84E  OO, 

2010  1  . S9303236E  OO,  .56651042E  OO,  .5393881  IE  OO,  .S1169418E  OO, 

2020  1  . 48345797E  OO,  . 45470942E  00,  .42547B99E  00,  . 39579765E  OO, 

2030  1  . 36S69686E  OO,  .335208S2E  00,  . 30436494E  00,  .27319881E  OO, 

2040  1  . 24174316E  00,  .21003131E  OO,  . 17809688E  OO,  . 1459737 IE  OO, 

2050  1  . U369S85E  OO,  . B1297495E-01 ,  . 48812985E-01 ,  . 16276745E-01  / 

2060  DATA  W/. 7967920 7E-03 , . 1B539608E-02 , . 291 073 1 8E-02 , . 39645S43E-02 , 

2070  1  . 50 1 42027E-02 ,  . 60S854S5E-02,  . 7096470BE-02 ,  .B1268769E-02, 

2080  1  .91 4867 1 2E-02 ,  . 10160771E-01 ,  . 1 1 162102E-01 ,  . 121S1603E-01 , 

2090  1  . 13128230E-01 ,  . 140909 42E-01 ,  . 1503872 lE-Ol ,  . 1 5970563E-0 l , 

2100  1  . 168854SOE-01 ,  . 17782502E-01 ,  . 18660680E-01 ,  . 195190B1E-01 , 

2110  1  . 20356797E-01 ,  . 21 172940E-01 ,  . 21966644E-01 ,  . 22737070E-O1 , 

2120  1  . 23483399E— 0 1 ,  . 24204842E-01 ,  . 24900633E-01 ,  . 25570036E-01 , 

2130  1  . 26212341E-01 ,  . 26826867E-01 ,  . 274 1 2963E-0 1 ,  . 2797000BE-0 1 , 

2140  1  . 2849741 lE-Ol,  . 2899461 4E-01 ,  . 29461 090E-01 ,  . 29896344E-0 1 , 

2150  1  . 302999 15E-0 l,  . 30671376E-01 ,  . 31010333E-01 ,  . 313I6426E-01 , 

2160  1  . 31589331E-01 ,  . 318287S9E-01 ,  . 32034456E-01 ,  . 32206205E-01 , 

2170  1  . 32343B23E-01 ,  . 32447 164E-01 ,  . 325161 19E-01 ,  . 325506 14E-01/ 

2180  NPOINT  -  48 

2190  PSI  -  (O. ,0. ) 

2200  G  -  U/2. /FLOAT <N) 

2210  DO  2  I  -  1,N 

2220  Z  -  (O. ,0. ) 

2230  DO  1  K  -  1, NPOINT 

2240  ARBA  -  G*(A(K)  ♦  2*1-1.) 

2250  ARGB  -  G* (-A (K)  ♦  2*1-1.) 

2260  Z  -  Z  ♦  (FCN (D,R, ARGA)  *  FCN <D ,R , ARGB) > *W (K> 

2270  1  CONTINUE 

2280  C  WRITE (6,990) I , K, FCN  <D,R , ARGA) , FCN (D , R , ARGB) 

2290  PSI  »  PSI  Z 

2300  2  CONTINUE 

2310  PSI1  -  G*PSI 

2320  990  FORMAT (2X,  ' I-  ,I2,2X,  K- ' , 12, 2X ,  F (A) -  * ,2E10. 4,2X , 'F (B>»* ,2E10. 

2330  PSI  -  <0. ,0.) 

2340  G  -  -8 

2350  DO  4  I  >  t,N 


Z  -  <0.  ,0.  ) 

DO  3  K  -  1, NPOINT 
ARGA  -  G*  < A  ( K )  2*1-1.) 

ARGB  -  G* (-A(K)  +  2*1-1.) 

Z  -  Z  +  < FCN (D,R, ARGA)  +  FCN (D,R, ARGB) >  *W(K) 

CONTINUE 

PSI  -  PSI  ♦  Z 

CONTINUE 

PS I 2  -  G*PSI 

PSI  -  PSI1  -  PSI2 

RETURN 

END 


2360 

2370 

2380 

2390 

2400 

2410  3 

2420 

2430  4 

2440 

2450 

2460 

2470 
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Z  -  <0.  ,0.  > 

DO  3  k  »  1,NP0INT 
ARGA  =  G* ( A (K)  +  2*1-1.) 

ARGB  =  G* ( —A  <K  >  *  2*1-1.) 

Z  =  Z  *  (FCN <D ,R, ARGA)  +  FCN < D ,R , ARGB) ) *W (K) 

CONTINUE 

PSI  =  PSI  ♦  Z 

CONTINUE 

PS  1 2  -*  G*PSI 

PSI  «  PSI1  -  PSI2 

RE  TURN 

END 
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LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  Investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  space  systems.  Versatility 
and  flexibility  have  been  developed  to  a  high  degree  by  the  laboratory  person¬ 
nel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  systems.  Expertise  in  the  latest  scientific  developments  Is 
vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The  labora¬ 
tories  that  contribute  to  this  research  are: 

Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  environmental  hazards,  trace  detection;  spacecraft  structural 
mechanics,  contamination,  thermal  and  structural  control;  high  temperature 
thermomechanics,  gas  kinetics  and  radiation;  cw  and  pulsed  laser  development 
including  chemical  kinetics,  spectroscopy,  optical  resonators,  beam  control, 
atmospheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmo¬ 
spheric  optics ,  light  scattering,  state-specific  chemical  reactions  and  radia¬ 
tion  transport  in  rocket  plumes,  applied  laser  spectroscopy,  laser  chemistry, 
laser  optoelectronics,  solar  cell  physics,  battery  electrochemistry,  space 
vacuum  and  radiation  effects  on  materials,  lubrication  and  surface  phenomena, 
thermionic  emission,  photosensitive  materials  and  detectors,  atomic  frequency 
standards,  and  environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  intelligence  and 
microelectronics  applications. 

Electronics  Research  Laboratory:  Microelectronics,  GaAs  low  noise  and 
power  devices,  semiconductor  lasers,  electromagnetic  and  optical  propagation 
phenomena,  quantum  electronics,  laser  communications,  Ildar,  and  electro- 
optics;  communication  sciences,  applied  electronics,  semiconductor  crystal  and 
device  physics,  radiometric  Imaging;  millimeter  wave,  microwave  technology, 
and  RP  systems  research. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metal 
matrix  composites,  polymers,  and  new  forms  of  carbon;  nondestructive  evalua¬ 
tion,  component  failure  analysis  and  reliability;  fracture  mechanics  and 
stress  corrosion;  analysis  and  evaluation  of  materials  at  cryogenic  and 
elevated  temperatures  as  well  as  In  space  and  enemy-induced  environments. 

Space  Sciences  Laboratory:  Magnet os pherlc,  auroral  and  cosmic  ray  phys¬ 
ics,  wave-particle  Interactions,  raagnetospherlc  plasma  waves;  atmospheric  and 
ionospheric  physics,  density  and  composition  of  the  upper  atmosphere,  remote 
sensing  using  atmospheric  radiation;  solar  physics.  Infrared  astronomy, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  Ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
Instrumentation. 


